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ABSTRACT 

We describe an object in the Hubble Deep Field North with very unusual near- 
infrared properties. It is readily visible in Hubble Space Telescope NICMOS images at 
1.6/im and from the ground at 2.2//m, but is undetected (with signal-to-noise < 2) in 
very deep WFPC2 and NICMOS data from 0.3 to l.l^m. The f^, flux density drops by 
a factor > 8.3 (97.7% confidence) from 1.6 to 1.1/im. The object is compact but may 
be slightly resolved in the NICMOS 1.6/Ltm image. In a low-resolution, near-infrared 
spectrogram, we find a possible emission line at 1.643^m, but a reobservation at 
higher spectral resolution failed to confirm the line, leaving its reality in doubt. We 
consider various hypotheses for the nature of this object. Its colors are unlike those of 
known galactic stars, except perhaps the most extreme carbon stars or Mira variables 
with thick circumstellar dust shells. It does not appear to be possible to explain 
its spectral energy distribution as that of a normal galaxy at any redshift without 
additional opacity from cither dust or intergalactic neutral hydrogen. The colors 
can be matched by those of a dusty galaxy at z > 2, by a maximally old elliptical 
galaxy at z > 3 (perhaps with some additional reddening), or by an object at z > 10 
whose optical and 1.1/xm light have been suppressed by the intergalactic medium. 
Under the latter hypothesis, if the luminosity results from stars and not an AGN, the 
object would resemble a classical, unobscured protogalaxy, with a star formation rate 
> IOOMq yr^-*^. Such UV-bright objects are evidently rare at 2 < z < 12.5, however, 
with a space density several hundred times lower than that of present-day L* galaxies. 

Subject headings: early universe — galaxies: photometry — galaxies: evolution — 
infrared: galax;ies 
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1. Introduction 

In recent years, astronomers have extensively developed the art of selecting interesting, high 
redshift objects on the basis of their broad band colors. Color selection has long been used 
to separate distant quasar candidates from the multitude of foreground stars either via their 
UV excess (e.g., Veron 1983) or by color criteria based on the passage of the 1216A Lyman a 
forest and 912A Lyman limit discontinuities through broad passbands (e.g.. Warren et al. 1987). 
Guhathakurta, Tyson, &; Majewski (1990) applied the latter technique to search for high redshift 
galaxies, and the method was brought to fruition by Steidel et al. (1996a,b, 1999) who have 
successfully identified and spectroscopically confirmed nearly 1000 galaxies at 2 < z < 4.5 using 
this approach. Other interesting classes of objects have been selected and studied on the basis 
of having extremely red optical-to-infrared colors (e.g., Elston, Rieke, &: Rieke 1988, 1989; Hu 
&; Ridgway 1994; Graham & Dey 1996; Thompson et al. 1999a). Some of these turn out to be 
early type galaxies at high redshift, while others are both distant and dust-obscured. In general, 
estimating galaxy redshifts from broad band colors is now a popular industry, and much progress 
has been made in applying such techniques at intermediate redshifts (e.g., Brunner, Connolly & 
Szalay 1999) and for identifying candidates for extremely distant galaxies (Lanzetta, Yahil, & 
Fernandez-Soto 1996). Indeed, some of the most distant galaxies now known have been found in 
this way (e.g., HDF 4-473, Weymann et al. 1998; HDF 3-951, Spinrad et al. 19980). 

The Hubble Deep Fields (North and South, or HDF-N and HDF-S, WilUams et al. 1996 
and 1999) have been valuable data sets for exploring such techniques because of the extremely 
deep, multiwavelength data which are available. The HDF-N has recently been the subject of two 
near-infrared surveys with NICMOS (Thompson et al. 1998) on board the Hubble Space Telescope 
(HST): a deep image of a ~ 1 arcmin^ sub-region (Thompson et al. 1999b), and a shallower 
map of the entire field (Dickinson et al. 1999; see also Dickinson 1999). These data offer new 
opportunities for identifying interesting objects on the basis of their colors. Here we describe an 
object with perhaps the most unusual colors in the HDF-N, which is significantly detected only 
at A > 1.6//m. The object was first noted by Lanzetta, Yahil, & Fernandez-Soto (1998), who 
used the ground-based infrared images which we obtained at Kitt Peak National Observatory (cf. 
Dickinson 1998) to identify candidate iT-band sources without optical counterparts. We discuss 
our imaging and spectroscopic observations of this object, and consider various interpretations of 
its nature. 

2. Imaging and Photometry 

We observed the HDF-N with NICMOS between UT 1998 June 13 and June 23, during the 
second refocus campaign when the HST secondary mirror was moved to ensure optimal focus for 



^We will occasionally use catalog numbers from Williams et al. 1996 to identify HDF galaxies in this paper. 
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NICMOS Camera 3. The observations and data reduction will be described in detail elsewhere 
(Dickinson et al. 1999); the immediately relevant aspects are summarized here. The complete 
HDF-N was mosaiced with 8 sub-fields, each imaged during three separate visits. During 
each visit, exposures were taken through both the FllOW (l.l^m) and F160W (1.6^m) filters. 
(Henceforth we will refer to the six WFPC2 and NICMOS HDF bandpasses as C/300, -B450, VeoGj 
-^8i4i <^iio and -f^ieo-) Each section of the mosaic was dithered through 9 independent positions, 
with a net exposure time of 12600s per filter, except in a few cases where telescope tracking was 
lost due to HST Fine Guidance Sensor failures. The region of interest for this paper did not lose 
any exposure time. The data were processed using a combination of STScI pipeline routines and 
custom software, and were combined into a single mosaic, accurately registered to the HDF-N 
WFPC2 images, using the "drizzhng" method of Fruchter & Hook (1999). The NICMOS images 
have a point spread function (PSF) with FWHM ~ 0^'22, primarily limited by the pixel scale 
(0''2) of Camera 3. Sensitivity varies over the field of view due to variations in NICMOS quantum 
efficiency and exposure time, but on average the images have a signal-to-noise ratio S/N ~ 10 
within an 0''7 diameter aperture at AB ~ 26.1 for both the Jno and Hiqq filters]^ In order 
to ensure properly matched photometry between the optical and infrared images, the WFPC2 
data were convolved to match the NICMOS PSF. Photometric catalogs were constructed using 
SExtractor (Bertin & Arnouts 1996), by detecting objects in the NICMOS images and measuring 
fluxes through matched apertures in all bands. 

We noticed the object which we will call HDF-N J123656.3+621322, or HDFN-JDl for 
brevity, during an initial visual inspection of the NICMOS data for objects with unusual colors. 
It is prominent at 1.6//m, but apparently invisible through all other HST filters, including Jno 
(Figure 1). Figure 2 shows a Jno — -f^ieo vs. -ffieo color-magnitude diagram highlighting the 
object's extreme and unusual color. In Jno — -H^ieoi HDFN-JDl is by far the reddest among the 
~ 1700 objects detected in our NICMOS survey. 

Careful inspection of the NICMOS data convinces us that this is a real astronomical source, 
not an artifact, and that there is no evidence that it is transient or variable. It is visible in each 
of the nine individual, dithered Hiqq exposures, which were taken during telescope visits on UT 
1998 June 16, 20 and 22-23, but is not detected in any of the corresponding Jno exposures taken 
during those same visits. We know of no NICMOS anomaly or optics ghost which could produce 
an artifact resembling what we observe.j] Moreover, the object is detected in two independent, 
ground-based near-IR data sets. It is faintly visible {S/N ~ 4) in the Kg (2.16/im) HDF images 
which we obtained with IRIM on the KPNO 4m telescope in 1996 April-May (cf. Dickinson 1998). 
HDFN-JDl is the brightest of the five -ftT-band selected, optically invisible candidates identified 



^Unless otherwise stated, we use AB magnitudes throughout this paper, defined as AB — 31.4 — 2.5 log(/,y), where 
(/^) is the flux density in njy averaged over the filter bandpass. 

^One NICMOS anomaly, the "Mr. Staypuft Effect" (cf. Skinner et al. 1998), produces electronic ghost images, 
but these are always offset by 128 pixels from brighter sources, which is not the case here. 
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by Lanzetta et al. (1998) from the KPNO IRIM HDF images. Their other fom candidates have no 
counterparts in the NICMOS Jno or Hiqq data, implying either that they are exceptionally red at 
A > 1.8//m, or that they are not real. 

To verify the Kg detection and to improve the quality of the photometry, we obtained new 
Kg images using NIRC (Matthews & Soifer 1994) on the Keck I telescope on UT 1999 April 5. 
108 dithered 60 second exposures were taken through light, intermittent cirrus and processed 
using the DIMSUM0 reduction software. The individual frames were scaled and weighted before 
combination using measured counts of a moderately bright, nearby star centered in the NIRC field. 
The combined image has a PSF FWHM = 0''5. The photometric zeropoint was bootstrapped from 
the IRIM Kg data using large aperture measurements of two stars, which gave excellent internal 
agreement. HDFN-JDl is readily visible as a compact source in the NIRC image (Figure 1). 

We summarize photometry of the object in Table 1. In order to set limits in the HST bands 
U300 through Jiio, we measured fluxes in a 0'.'7 diameter circular aperture at the -ffieo centroid 
position. This aperture maximizes S/N for the object in the -ffieo image, and is very close to that 
required to maximize S/N for point sources. We also have verified that the optical measurements 
in Table 1 are consistent with limits derived from the original, unconvolved WFPC2 images. A 
21% aperture correction, derived from the i^ieo image, has been applied to the measurements, 
errors and limits in all HST bands to adjust the fluxes to "total" values. This correction assumes 
that the object has similar morphology at all wavelengths, which cannot be verified at present. 
We followed a similar procedure for the NIRC Kg image, using photometry measured within a 1''2 
diameter aperture, corrected to total flux using curve of growth measurements of the reference 
star. We also measured the object in the IRIM Kg image using a method based on that of 
Fernandez-Soto, Lanzetta & Yahil (1999). The NICMOS image is used as a template, convolved to 
match the IRIM PSF, and then fit to the ground-based data. The NIRC and IRIM measurements 
{Kab = 23.9 and 23.8, respectively) agree within their errors, and with the value measured by 
Lanzetta et al. (1998) (Kab = 23.7) from the same IRIM data. The la and 2a color limits for 
Jiio — HiQo are >2.8 and >2.3, respectively, while -ffieo — Kg = 1.2'i~^^'^. 

HDFN-JDl has S/N < 2 in all bands other than Hiqq and Kg, regardless of the aperture size 
used. However, we do measure formally positive counts above the background at -B450, Vqqq, I^u 
and Jiio- Although measurements at low significance levels are prone to systematic uncertainties 
in background subtraction, etc., there is nevertheless the possibility that the object does have 
non-zero optical flux. To explore this further, we have used an additional 63000s WFPC2 I814 
image of the HDF obtained by Gilliland, Nugent, & Phillips (1999) in a search for high redshift 
supernovae. R. Gilliland kindly provided his sum of all available HDF-N I814 data in the WF3 
region, with a total exposure time of 186600s. We registered this image with the data in other 
bandpasses. Photometry on the unconvolved, "grand sum" /814 image in an 0''7 aperture centered 



*^Deep Infrared Mosaicing Software, a package of IRAF scripts by Eisenhardt, Dickinson, Stanford and Ward, 
available at ftp://iraf.noao.edu/contrib/dimsumV2. 
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at the nominal position of JDl measures a flux of 5.7 it 2.6 nJy (with no aperture correction), or a 
formal S/N = 2.2. 

To assess the significance of this optical measurement, we carried out a simple fluctuation 
analysis on the WFPC2 data. The images were normalized to constant variance over the fleld 
of view by dividing by smoothed noise maps which are generated as part of the data reduction 
process. We used SExtractor to generate masks which exclude all readily detectable sources in 
each image, including generous buffer regions around each. We fit smooth background maps 
to the "blank sky" regions and subtracted them from the images. The data were filtered by a 
Gaussian with FWHM = 0''14 (matching the WFPC2 PSF), and the distribution of pixel values 
over an 855 arcsec^ region was compared to that near the position of HDFN-JDl. We did this 
for each WFPC2 image, for variance weighted sums of the optical bandpasses, and also for 
combinations constructed following the procedure of Szalay, Connolly, & Szokoly (1999). The 
results from the images were statistically similar to those from the weighted sums, so we 
describe only the latter here. The filtered pixel values have a nearly Gaussian distribution (whose 
width in a given image we will characterize as a here), with a positive tail due to faint sources. 
At the nominal JDl peak position, the filtered /814 image has a value which exceeds the local 
background by l.Gcr, and in the weighted Vqqq + Isu image by 1.9a. There is a stronger local 
peak in the filtered I814 and V606+-^8i4 images located 0"14 from the nominal JDl position, which 
exceeds the local background by 3.1a in /8i4 and by 3.3a in Vqqq + 1814. Including the -B450 data 
in the weighted sums or images does not strengthen the peak. Its formal significance depends 
on how one treats the non-Gaussian tail of positive pixel values due to faint sources. Excluding 
buffered regions around all objects detected in the infrared catalogs, we find that there is a 5.1% 
chance of a pixel in the filtered V606+-^8i4 images exceeding this "3.3o"" threshold within O'.'lA of 
a location specified a priori. The infrared catalogs, however, miss faint blue galaxies which are 
readily visible in the optical data, and which contribute to this positive tail. Excluding all objects 
detected in the infrared or optical SExtractor catalogs, the chance probability drops to 1.8%. The 
exact probabilities depend on the detection threshold used for the catalogs. We conclude that the 
optical "detection" corresponds to a Gaussian significance of 1.6 to 2.2cr. 

HDFN-JDl is not detected in the VLA radio maps of Richards et al. (1998) and Richards 
(1999), with a 3a limit of 4.8/xJy at 8.5 GHz and 22.5/xJy at 1.4 GHz. It is not reported as a source 
in the ISO HDF catalogs of Goldschmidt et al. (1997), Aussel et al. (1999) or Desert et al. (1999), 
with approximate flux limits of 50;uJy and 25/xJy at 7 and 15/xmJ^ nor as an 850/xm source in the 
SCUBA observations of Hughes et al. (1998), with a limit of 2.0 mJy (4.4o"). Interestingly, the 
next reddest HDF object, J123651. 74+621221.4 (with Jno - Hieo = 1.6, i^ieo = 24.3; see Figure 
2), corresponds to a VLA/MERLIN radio source, and may also have a 15/im counterpart in the 
Aussel et al. (1999) ISO source list, as well as a possible 3a detection in the 1.3mm IRAM map of 



^These ISO flux limits are estimated based on measurements reported by Aussel et al. (1999) for two other 
"supplemental" 15/xm sources in the same general region of the HDF. 
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Downes et al. (1999). 

3. Spectroscopy 

We observed HDFN-JDl using the Cryogenic Spectrograph (CRSP; Joyce 1995) at the 
KPNO 4m on UT 1999 May 2, 3, 6 and 7. We acquired the target using a bhnd offset from a 
bright star located in the HDF flanking fields. The observations were taken through a l!'0 wide 
slit oriented at PA = 119°4 to cover both the target object and a nearby star (12''8 away) which 
served as a pointing reference. The target was dithered along the slit in an ABBA pattern using 
100 second exposure times, and the count rate from the reference star was monitored to ensure 
that pointing remained stable. The target was periodically reacquired and placed at new positions 
along the slit. Seeing throughout all observations averaged l!'0. On 2 and 6 May we used grating 
4 (200 lines/mm, blazed at 3;um) at 1st order in the i^-band (AA1.90-2.50^m) giving an effective 
spectral resolution R = A/AA 240 at A2.2/im as measured from night sky lines. On 3 and 6 
May we also observed with grating 4 at 2nd order in the iJ-band (AAl.49-l.80/xm), with R k, 380 
at A1.65/Ltm. The iJ-band data on 3 May were taken through occasional cirrus; frames with poor 
transparency (judged from the reference star) or bad sky subtraction were discarded from the 
analysis. The total exposure times retained for the final H and i^-band grating 4 spectrograms 
were 13400s and 17600s, respectively. We reduced the data using standard procedures. The data 
were corrected for array non-linearity, and flat-fielded using dome flats. The sky was subtracted 
(to first order) using the ABBA differences. The wavelength scale and geometric distortion were 
calibrated using OH night sky lines, and the two dimensional spectral images were rectified. 
Residual sky features were then fit and subtracted. The images were positionally registered using 
the reference star, and combined using a bad pixel mask and a percentile-clipping scheme to reject 
outlying pixels. Spectra were extracted through a 3 pixel (l''05) wide window at the nominal 
position of the object relative to the reference star, along with a noise measurement extracted 
from the 2-dimensional variance image created when the images were combined. Flux calibration 
was based on observations of the Elias et al. (1982) standard HD 105601, with absolute flux 
normalization based on the H and Ks magnitudes (measured from the IRIM HDF data) of the 
reference star on the slit. 

The low-resolution iJ-band spectrogram (Figure 3) shows a possible emission line at 
A1.643/im, with flux ^ 2 y. 10~^^ erg s~^ cm~^ and integrated S/N ~ 4.3. Although the i^-band 
sky is nearly covered with OH emission bands at this low dispersion, the telluric lines near this 
wavelength are relatively weak. We reobserved the object at higher spectral resolution (i? k, 670 at 
A1.65/xm) on UT 1999 May 7 under good atmospheric conditions, using grating 1 (300 lines/mm, 
blazed at Ajim) in 2nd order, covering AA1.51-1.69jLtm. The putative emission line does not 
reproduce in the grating 1 spectrogram. The 3cr upper limit for an unresolved emission line at 
1.643/xm is f« 6 X 10~^^ erg s~^ cm~^. 

Careful, frame-by-frame inspection of the grating 4 data does not reveal any obvious artifacts 
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that could have produced the emission feature. It is present when the data are combined as a 

median (rather than percentile clipped averages), demonstrating that the feature results from 
systematically high data values and not from intermittent artifacts. We have tried to test its 
reality by subdividing the grating 4 exposures into independently averaged, randomly selected half 
data sets. However the S/N in each half-set is low; the apparent line sometimes appears in both, 
sometimes not. Its significance is hard to evaluate because noise in the data is correlated (e.g., sky 
subtraction residuals) with greatly varying amplitude vs. wavelength. As a test, we normalized the 
coaddcd, two dimensional spectral image by the variance map to equalize the pixel-to-pixel RMS 
throughout, and extracted 26 spatially independent, 3 pixel wide spectra from regions of the slit 
unaffected by the reference star. Each extraction was subtracted to zero mean, and then convolved 
by a Gaussian with the instrumental resolution. We measured the peak value of the nominal 
emission line from JDl, and searched for features with equal or greater amplitude (positive or 
negative) at any wavelength in the other 25 "test" extractions. No comparable positive fcatTirc 
was found, and only one negative feature, suggesting a probability ~ 2% (1 out of 2 x 25) that the 
HDFN-JDl "line" would arise by chance. 

Our spectroscopic results are therefore ambiguous. The line in the grating 4 spectrogram 
is resilient but not ironclad. The fact that it does not reproduce in the grating 1 observation 
suggests either that it is not real, or that it is well resolved at the higher dispersion (thus reducing 
the detection sensitivity), or that the object was not well placed in the spectrograph slit on May 7. 

The -fC-band spectrogram shows no significant emission features; the 3a flux limit for an 
unresolved line ranges from 1 to 3 x 10~^^ erg s""^ cm~^ over the range AA1.95-2.4//m, and rises 
steeply at longer wavelengths due to the increasing thermal background. We would not have 
expected to (and did not) detect the object's continuum at either K or H. 

4. Angular Extent 

The NICMOS F160W image of HDFN-JDl is very compact, but there is evidence that it 
may be resolved. Examining ten well exposed, spectroscopically confirmed stars in our images, 
we measure the PSF FWHM to be 0''217 ± 0''016 (0''19 to 0''24 maximum range), while JDl has 
FWHM = 0('28. Figure 4a compares the surface brightness profiles of the object and three faint, 
well isolated stars which have been registered and scaled using a non-negative cross-correlation 
procedure. The profile of HDFN-JDl appears to be slightly more extended. Subtracting the 
scaled PSF stars from JDl leaves positive annular residuals which are not present when one star 
is scaled and subtracted from another. Additionally, SExtractor computes a stellarity index using 
a neural network classifier which outputs values from (extended) to 1 (stellar). The known HDF 
stars have stellarity > 0.84, while JDl has stellarity = 0.08 (Figure 4c). However, the object is 
near the faint limit where the classifier appears to be reliable. The NICMOS PSF may depend 
on source color, and HDFN-JDl is extremely red while the stars are among the bluest objects 
in Jiio — HiQQ. We have examined this using Tiny Tim models (Krist &; Hook 1997) computed 
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for a very wide range of source spectra and find it is too weak to account for the differences we 
measure. Proper analysis of the PSF from dithered images with NICMOS Camera 3 requires a 
reahstic treatment of the known (and large) detector intrapixel sensitivity variations (cf. Lauer 
1999). Therefore while our measurements suggest that the object is spatially extended, we cannot 
be completely confident about this given the present data. In the future, a robust angular extent 
measurement might be accomplished using NICMOS Camera 2, which critically samples the PSF 
at 1.6/im, when the instrument is resuscitated with the NICMOS Cooling System in HST Cycle 10. 

Figure 5 compares the colors of HDFN-JDl to those of a variety of cool or reddened stars. 
We have converted our NICMOS AB magnitudes to a conventional Vega scale for comparison 
to published stellar data, including an approximate color term to correct to standard J and H 
bandpasses. In this system, JDl has J — H > 2.5 {2a) and H — K 1.6. This is much redder than 
ordinary cool stars and known substellar objects, which can have extremely red optical-to-infrared 
colors but are generally rather blue in J — H and H — K where their spectra are dominated by 
strong molecular absorption. The reddest L dwarf stars reach J — H = 1.45, but the known 
"methane dwarfs" have J — H^H — K<0. Atmosphere models (Burrows et al. 1997) predict 
still bluer colors for cooler brown dwarfs and giant extra-solar planets. The only known stars as 
red as JDl are those undergoing mass loss with thick circumstellar dust shells, such as extreme 
carbon stars or Mira variables. HDFN-JDl is redder than ordinary galactic or Magellanic AGB 
stars, but the most heavily reddened objects like IRC +10216 and some extreme Miras can 
equal or exceed its colors .|^ It would be remarkable to encounter such an unusual star in a tiny, 
high galactic latitude field such as the HDF. Moreover, at K it is 11.9 magnitudes fainter than 
the reddest LMC carbon stars, implying (by analogy) a distance ~ 12 Mpc. A blackbody with 
T ~ 1050K also matches the -f^ieo ~ Kg color and 2ct Jhq — -ffieo limit for HDFN-JDl quite well, 
but no known star has such a spectrum. Recent atmosphere models for cool white dwarfs (e.g., 
Hansen 1998; Saumon & Jacobson 1999) suggest that they should have very blue colors due to H2 
opacity. 

5. Discussion 

The salient photometric features HDFN-JDl are firm detections at i?i60 and Ks, and very 
faint upper limits (possibly with some S/N < 2 detections) at all shorter wavelengths. The 
bandpass-averaged fy flux density declines by a factor of 3.1 from 2.2 fim to 1.6/im, and then again 
from 1.6/im to l.l//m by a factor 34^^ {la errors), with a 97.7% (1-sided) confidence limit > 8.3, 
implying curvature or a break in the spectral energy distribution. This is the only object in the 
NICMOS HDF-N at Higq < 26 which is undetected (with S/N < 2) in any of the optical WFPC2 
bandpasses. It is by far the reddest HDF-N object in Jno — Hiqq, but also one of the reddest 



IRC +10216 itself, scaled to the -ffieo or Ks magnitude of HDFN-JDl, should have been detected in the ISO 
images of the HDF at 7 and 15/^m. 
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in Hi6o — Kg. Only one (slightly) brighter HDF object, a z = 3.2 "i7-dropout" galaxy, has a 
redder .ffieo — I^s color (but equal within the la measurement uncertainties). That galaxy's color 
probably results from a strong Balmer and/or 4000 A break, or possibly from strong [0III]+H/3 
line emission in the Kg band. 

There arc three common reasons why a galaxy's spectral energy distribution (SED) may 
appear to be very red: extinction, age, or the presence of a strong spectral break, such as that 
caused by the Lyman limit or Lyman a forest. Here we consider each of these in turn, making 
comparisons wherever possible to photometry of actual galaxies from the HDF and elsewhere 
in order to avoid overreliance on models. One way of comparing colors in different bandpass 
combinations is to parameterize the photometry by a spectral index a, i.e. fu oc For two 

bandpasses with effective wavelengths Ai and A2 and magnitude measurements mi and r?i2, we 
may define CK(mi— ~ 

0.4(mi - m2)/log(A2/Ai). The nominal HDFN-JDl Jno - i?i60 color 
corresponds to a(^j_H) = 10.2, and the la and 2a limits to q:(j_h) > 7.5 and > 6.1, respectively. 
Considering pairs of adjacent HDF bandpasses [/300 through Jno, we find a few other galaxies with 
comparably steep spectra at shorter wavelengths. Among these, the objects with spectroscopic 

identifications are invariably either Lyman break "dropouts," or red, early-type galaxies, and the 
as-yet unidentified objects all appear to be consistent with being members of one of these two 
classes. 

Dust may, in principle, redden a spectrum almost arbitrarily. Because extinction laws 
generally steepen in the UV, a red and possibly strongly curved SED may result when the UV 
portion of a distant, dusty galaxy redshifts into optical or near-IR bandpasses. One well-known 
and fairly extreme example is the object discovered by Hu Sz Ridgway (1994) and colloquially 
known as HRIO, with z = 1.44 (Graham & Dey 1996). Sub-mm measurements by Cimatti et 
al. (1998) and Dey et al. (1999) have shown that it is probably a dust-enshrouded, star forming 
galaxy or AGN. The SED of HRIO is steepest between / and J (Dey et al. 1999), with a(j_j) = 6.1. 
An object like HRIO shifted to z ~ 2.3 would have colors roughly consistent with our measured 2a 
limits on HDFN-JDl (sec Figure 6). A somewhat redder SED might be needed to fully match the 
photometry of JDl, but this could presumably be accomplished by adding still more dust. HRIO is 
compact in Keck iC-band images, but diffuse (f« l''2 in size) and bimodal in WFPC2 images which 
sample the rest-frame UV light (Dey et al. 1999). Similar changes in morphology with wavelength 
could make it difficult to detect a fainter, more distant analog to HRIO in the WFPC2 HDF. 
Adopting (here and below) a cosmology with Q.m = 0.3, Oa = 0.7, and Hq = 70 km s~^ Mpc~^, at 
z 2.3 HDFN-JDl would have a rest-frame 5-band luminosity 0.12x that of HRIO. Scaling 
the far-IR emission from HRIO accordingly, the predicted 850/xm flux would be ^0.5 mJy, or a 
factor of four below the current SCUBA limits for the HDF. Objects as red as HRIO are relatively 
rare, and it may seem surprising to find an even redder example within the 5 arcmin^ of the HDF, 
but little is known about their numbers at these faint magnitudes. The fact that there are only a 
handful of objects with such extreme IR colors in the HDF suggests that they are not a common 
population by number relative to ordinary, relatively unreddened galaxies, but their unobscured 
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luminosities and star formation rates might be quite large and important in the scheme of galaxy 

formation. Some of the faint sub-mm sources detected in recent surveys arc red objects like HRIO 
(e.g., Smail et al. 1999), although it is worth noting that none of the candidate counterparts to 
sub~mm sources in the HDF has particularly unusual colors (Dickinson et al. 1999). 

The integrated spectrum of an old stellar population is steepest in the near ultraviolet, from 
the famihar 4000A break through spectral breaks at 2900A, and 2640A (Morton et al. 1977). 
At 2.5 < z < 4.5 these breaks would redshift beyond the WFPC2 bandpasses into the near-IR, 
maximizing the Jno — i^ieo color. The resulting /c-correction could make old, high redshift 
ellipticals nearly invisible in the optical HDF (e.g., Maoz 1997). For example, the observed-frame 
colors of the giant elliptical galaxy HDF 4-752.1 at z = 1.013 correspond to c((b-v) ~ 10-2 and 
^{V-i) ~ 6.3. Shifted to z > 3, its SED would roughly match that of HDFN-JDl (Figure 6). 
However, this would require that a z « 3 elliptical have a spectrum nearly as red as that of 4-752.1 
at z = 1, despite the universe being considerably younger. This is difficult to explain without 
invoking an unusual IMF, an unpopular cosmology, or extinction. Figure 7 shows the infrared 
colors of HDF galaxies with known redshifts, along with models computed using the population 
synthesis code of Bruzual & Chariot (1993, 1996). The reddest model traces a solar metallicity 
population formed with a Salpeter IMF in a short burst at z = 15 and aging passively thereafter. 
The HiQo — Kg color of HDFN-JDl is reasonably matched for 3 < z < 6, but the Jno — -ffieo color 
of the model does not quite reach the formal 2a color limit for JDl at any redshift. At z = 3.5 
in our adopted cosmology, JDl would have an absolute B magnitude Mb ~ — 21.6. If it were an 
"old" elliptical galaxy at z > 3, this would strongly suggest that at least some galaxies formed 
the bulk of their stars at extremely large redshifts. If so, however, then it is puzzling that there 
are no other HDF galaxies nearly as red, as might be expected if there were a continuum of such 
objects extending out to z ^ 3. Red, early-type galaxies have been identified spectroscopically 
out to z = 1.55 (e.g., Dunlop et al. 1996; Spinrad et al. 1997; Dickinson 1997; Dey 1998) and 
photometrically in the HDF and elsewhere out to 2; 2 (cf. Stiavelli et al. 1999, Benitez et 
al. 1999), but HDFN-JDl would stand alone as a unique example at significantly larger redshift. 

The Lyman limit can effectively truncate the spectrum of high redshift objects, making them 
appear arbitrarily red in certain color combinations. The z > 5 "V— dropout" galaxies HDF 
4-473 and 3-951 (Weymann et al. 1998, Spinrad et al. 1998) and the z = 4.022 "B-dropout" 
HDF 3-512 (Dickinson 1998) all have a > 7.3 for colors that span the Lyman limit, and several 
2; ~ 3 "U-dropouts" have a^jj-B) > 6. At z > 3, an additional break is introduced by the 
Lyman a forest, which is increasingly thick at higher redshifts and affects broad band colors to a 
correspondingly greater degree (cf. Lowenthal et al. 1997, Dickinson 1998). We do not know the 
opacity of the Lyman a forest at 2; > 6, but may reasonably assume that it continues to increase 
since we believe that the universe was dominated by neutral hydrogen beyond some reionization 
redshift. 

If we interpret HDFN-JDl as a Lyman break object, then its Jno — i^ieo color is matched 
for z > 10, but the red -ffieo ~ Ks color suggests that the Lyman a forest has eaten away roughly 
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half of the i?i6o Aux, implying a redshift z ^ 12.5 (see Figure 7). At that redshift, the K-band 

corresponds to Aq = I6OOA in the emitted frame, essentially the same wavelength where Steidel 
et al. select galaxies in their survey at z ~ 3. We may therefore make a direct comparison to 
the observed 2: « 3 luminosity function (Dickinson 1998, Steidel et al. 1999). For an Q^m = 0.3, 
Oa = 0.7 cosmology, HDFN-JDl would be 3x more luminous than the brightest Lyman break 
galaxy from our 0.3 deg^ ground-based survey. This may seem improbable, given the small solid 
angle of the HDF. However, the total co-moving HDF volume out to z = 12.5 is quite large, 
particularly for an open or A-dominatcd cosmology. For Q,m = 0.3, = 0.7, it is equivalent to 
the effective volume of the UnGTZ Lyman break selection function (see Steidel et al. 1999) over a 
100 arcmin^ field. Because the volume is heavily weighted toward the highest redshifts, then it is 
not unlikely that the rarest, most luminous objects would also be the most distant, provided that 
they exist at all beyond z > 6. In our ground-based survey we often find 2; pa 3 QSOs that are 
brighter than the most luminous Lyman break galaxies, and the nearly unresolved morphology of 
JDl might plausibly indicate that it is some sort of AGN. Populations of distant QSOs have been 
postulated as a means of reionizing the universe at high redshift, although existing models for the 
formation of high redshift AGN (Haiman & Loeb 1998; Haehnelt, Natarajan, & Rees 1998) would 
predict that no objects as bright as JDl should be found within the HDF volume at 2; > 10. 

If HDFN-JDl were a galaxy forming stars at z = 12.5 with a Salpeter IMF and without dust, 
its UV luminosity would correspond to a star formation rate ISOh^^MQ/yi (70 or 4OO/ifo^M0/yr 
for Einstein-de Sitter or = 0.2 open universes). This is the sort of rate required by monolithic 
collapse models which would form a W^^Mq galaxy within a short (< lO'^ yr) time scale, and is 
also comparable to the obscured star formation rates which have been claimed for high redshift 
sub-mm sources. However, if HDFN-JDl resembles a classical, unobscured protogalaxy, then such 
objects are evidently quite rare. We will present a more complete discussion of color-selected high 
redshift galaxy candidates from our NICMOS survey in a future paper, but there are no U, B, V 
or I "dropout" candidates in the HDF with luminosities comparable to that which JDl would have 
if it were at z = 12.5. One such object in the redshift range 2 < z < 12.5 implies a space density 
~ 10~^hjQ Mpc^'^ for our assumed cosmology, and 4.3 x larger for an Einstcin-dc Sitter universe. 
This is a few hundred times smaller than the present day space density of L* galaxies {(f)* = 4 
to 6/170 ^ 10"^ Mpc^'^ from the i^-band luminosity functions of Mobasher et al. 1993, Gardner et 
al. 1997, Szokoly et al. 1998, and Loveday 1999). If rapid, monolithic galaxy formation took place 
anywhere in that redshift range, then either it was quite uncommon, or it was obscured by dust. 
Lanzetta et al. (1998) reached very similar conclusions based on their search for optically invisible 
objects in the KPNO infrared HDF data. If perhaps only the rarest, most massive galaxies, 
e.g., brightest cluster ellipticals (BCEs), formed their stars rapidly at the highest redshifts, then 
HDFN-JDl could be one example. The present-day space density of galaxy clusters with Abell 
richness class > 1 or X-ray temperature kT > 2.5 keV is 3 x IO^^^^q Mpc~^ (cf. Eke, Cole, & 
Frenk 1996), implying that we might expect ~ 1/3 "proto-BCEs" per HDF volume. 

From Figures 7 and 1 it is evident that a star-forming galaxy at 2; > 9.5 should have 
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Jiio — ^^160 ~ 1-5 due to the Lyman limit, and also that HDFN-JDl is the only such object in 
the HDF-N with Hiqq < 26, the magnitude down to which we could measure such a color or 
limit with > 2a significance.]^ There are, therefore, no other detected candidates for galaxies at 
z > 9.5 in the HDF, although cosmological surface brightness dimming could significantly affect 
sensitivity to extended protogalaxies at such redshifts, even in deep images such as these. At 
z = 9.5, the -ffieo = 26 limit corresponds to an unobscured star formation rate ^ IOHjqMq yr~^ 
for our adopted cosmology. 

The possible 1.643;^m emission line seen in the CRSP spectra could be consistent with any 
of these interpretations. For a dusty galaxy or AGN, it might be [OIIIJA5007A at z = 2.28. Ha 
at z ~ 1.5 is also possible, but it is much easier to match the colors at z > 2 where the Jno 
band would sample the rest-frame UV. Alternatively the line could be [OII]A3727A at z = 3.40, 
a redshift where an old stellar population nearly matches the colors of HDFN-JDl. An emission 
line, however, might suggest active star formation inconsistent with a maximally old elliptical 
galaxy, thus requiring either the presence of an AGN, dust, or both. A reddened, star forming 
object at this redshift is also possible. Finally, the line could be Lyman a at z = 12.51, where the 
JHK colors are well matched by the Lyman break hypothesis. However, given the fact that the 
line did not reproduce in our reobservation (see §3 above), we cannot be confident that it was real 
and do not further consider the spectroscopic possibilities here. 

An alternative explanation for the peculiar colors would be the presence of very strong line 
emission in one or more infrared bands. Line fluxes ~ 1.5 and 2 x 10~^^ ergs~^ cm~^ in the -ffi60 
and Ks bands, respectively, could mimic the broad band fiuxes measured for HDFN-JDl. The 
possible 1.643^m emission line, if real, could therefore account for the entire signal detected in the 
NICMOS -ffieo image. Detection at both i/ieo ^-nd Kg makes the "pure emission line" hypothesis 
seem less plausible, however, requiring either strong lines in both bands or an extremely red 
continuum. In the latter case, we are forced back to our previous speculations. 

6. Conclusion 

Without further data, we cannot unambiguously distinguish between these explanations for 
the nature of HDFN-JDl. Each, however, is quite remarkable in its own right. Perhaps the 
least "exotic" extragalactic hypothesis is that JDl is a dusty, HRlO-like galaxy at z > 2. Such 
objects, with sub-mm fluxes just below the current SCUBA detection limits, might be sufficiently 
common to make up the bulk of the far-infrared background (Barger, Cowie & Sanders 1999). 
The possibility that HDFN-JDl is an "old" elliptical galaxy at z > 3 is perhaps more remarkable, 
in that it would strongly suggest that at least some galaxies formed the bulk of their stars at very 



^^The radio source counterpart J123651. 7+621221. 4 is well detected in the -B450, Veoe and 78i4 optical passbands, 
and thus is not at z 2> 3. 
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large redshifts.0 

Perhaps the most spectacular interpretation of this object, but also the one requiring the 
most rigorous proof, would be that it lies at z > 10, an extremely distant analog of the Lyman 
break galaxy population. If this were the case, then its high luminosity and the relatively small 
volume of the HDF would suggest that we were either extremely lucky, or that such objects were 
more common than would be expected based on an extrapolation of the z ~ 3 population to 
higher redshifts. If the luminosity results from star formation, then HDFN-JDl would resemble 
the classical picture of a protogalaxy, forming stars at > lOOM0/yr. However, the implied space 
density of unobscured objects with such star formation rates is far smaller than that of L* galaxies 
today. If galaxies formed monolithically at high redshift, then this was either a rare occurrence, 
or the process was largely enshrouded by dust, as has been suggested by the recent detection 
of distant sub-mm sources. If, instead, JDl were a QSO, it could be part of the population of 
objects responsible for re-ionizing the universe. From our NICMOS HDF data, we should be able 
to identify unreddened, compact objects with star formation rates ~ 20h^f^MQ/yi at z ~ 9.5, with 
the limit rising to ~ lOOh^Q Mq/yt at z = 13 as the Lyman a forest suppresses the -ffieo flux. No 
candidates other than HDFN-JDl are seen. 

Distinguishing between the possible explanations for this object will require new observations, 
but none will be easy. Spectroscopy will be challenging even with the new generation of IR 
instruments on 8 to 10m telescopes, and there is no guarantee that the object has detectable 
emission lines. Although our understanding of the transparency of the IGM at z > 5 is based solely 
on extrapolation, a robust optical detection would probably exclude the Lyman break hypothesis. 
Occasionally, QSO sightlines at lower redshifts are free of optically thick H I absorption (e.g., 
Reimers et al. 1992), but it seems unlikely that emitted-frame Lyman continuum from an ordinary 
galaxy or QSO could propagate through the universe from z ~ 12 without being absorbed. The 
optical measurement in §2 is suggestive but not highly significant; a deeper optical image, e.g., 
with STIS or the forthcoming HST Advanced Camera for Surveys, would be a valuable (albeit 
expensive) observation. Photometry at A > 3/um with IRAC on SIRTF could distinguish between 
the z > 10 Lyman break and z < 6 red galaxy hypotheses, since the spectra of red galaxies should 
continue to rise toward longer wavelengths, with fu > 2/iJy, while a star forming galaxy at z > 10 
should have a flatter, fainter SED (see Figure 6). If near-IR spectroscopy fails to detect emission 
lines, then the best hopes for distinguishing between the old galaxy and reddened starburst/AGN 
models lie at longer wavelengths. Given that the HDF already has the deepest (and possibly 
confusion limited) SCUBA observation, a true sub~mm detection may require a future generation 
of telescopes and instruments. A detection would certainly imply dust, but given the negative 
sub-mm /c -correction it probably would not distinguish between "low" (z ~ 2) and high (z 12) 
redshifts without multifrequency measurements. For a reddened starburst at z ~ 2.3, the dust 
emission should peak near ~ 200^m and might be detectable by the SIRTF MIPS instrument 



When first forming stars, the 2: > 10 progenitor to an old, z « 3.5 elliptical might well resemble HDFN-JDl. 
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near its confusion limit. SIRTF IRAC photometry from 3.6 to 8fxm may prove useful, although 
disentangling reddening and age from broad band SEDs is notoriously difficult. 

Regardless of its nature, the fact that such an unusual and extreme object was found in 
a 5 arcmin^ field suggests that there are interesting surprises awaiting future, wide-field, deep 
infrared surveys. 
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Table 1. Photometry of HDF-N J123656.3+621322 



Instrument 


Bandpass 


Wavelength 


Flux Density'' 


AB Magnitude'' 


NIRC 




2.16/im 


995.0 ± 182.5 nJy 


23.911°:?^ 


NICMOS 




1.61/im 


318.1 ± 13.8 nJy 


or 14+0.05 


NICMOS 


JllO 


1.14/im 


9.3 ± 14.6 nJy 


> 27.44 


WFPC2 


-^814 


0.80/im 


6.1 ±4.1 nJy 


> 28.51 


WFPC2 


Vfioe 


o.eo^m 


4.3 ± 2.4 nJy 


> 29.00 


WFPC2 


^450 


0.46/jm 


1.5 ± 3.6 nJy 


> 29.05 


WFPC2 


t^300 


0.30/im 


-2.5 ± 7.5 nJy 


> 28.66 



''Fluxes and uncertainties with aperture corrections to "total" values (see text). 
''When S/N < 2, magnitudes are quoted as 2cr limits. 
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Fig. 1. — HST and Keck images of HDFN-JDl at 0.3-2. 16|Ltm. The field of view of each panel is 
4" X 8", and north is 23° 8 counterclockwise from vertical. HDFN-JDl is identified by tick marks 

in the Hiqq and Ks panels, and is located at J2000 coordinates a = 12:36:56.32, 6 = 62:13:21.7. 
The NIRC Kg image has been smoothed by a Gaussian with FWHM = 0''38. A constant grey level 
corresponds to constant /j/ surface brightness. The elliptical galaxy at top (HDF 3-48.0) is very 
red in its own right; its estimated photometric redshift is 2; 1.25 (Budavari et al. 1999). 
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Fig. 2. — Jiio — ^^160 versus i^ieo color-magnitude diagram (on the AB system) for the 
HDF-N NICMOS-selected catalog. Filled and open squares indicate galaxies with and without 
spectroscopic redshifts, respectively, while known stars are marked by star symbols. Colors are 
plotted with ±lcr error bars, and measurements with S/N < 2 are marked by triangles at the 2a 
color limit. HDFN-JDl is labeled; its la color limit would be Jno — -ffieo > 2.8. The next-reddest 
galaxy, J123651. 7+621221.4, is a faint cm radio source. Another VLA source, J123642.1+621331.5, 
is also unusually red (Richards et al. 1998, Waddington et al. 1999). 
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Fig. 3. — -ff-band spectrograms of HDFN-JDl obtained with CRSP at the KPNO 4m. The spectra 
have been smoothed by a 3-pixel boxcar. The hghtly shaded regions show the it la noise level of 
the data (before smoothing), and greatly suppressed sky spectra are plotted along the bottom axis. 
The low-resolution grating 4 spectrogram (top) shows a possible emission line at A1.643;um, but 
this is not detected in the subsequent, higher resolution grating 1 spectrogram (bottom), leaving 
its reality in doubt. 
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Fig. 4. — Left: (a) Radial surface brightness profile of HDFN-JDl (points) compared to three faint, 
isolated HDF stars (lines, shaded region) which have been registered and scaled by non-negative 
cross-correlation. Right top: (b) Hiqq half light radius vs. magnitude for HDF objects. Known 
stars are indicated by star symbols, and HDFN-JDl is labeled. Right bottom: (c) SExtractor stellar 
classifier (at -ffieo) for HDF objects, with symbols as in (b). 



- 25 - 



S) 2 

to 

> 



I 



T 



T p r 



T 



T 



T 1 r- 



-1 



+ LMC C+M AGB stars 
X IRAS Miras + Carbon stars 
A LMC red supergiants 
□ Bulge M giants 
M dwarfs 
L dwarfs 
Methane dwarfs 



IRC +10216 



X 




LHS 3250 



_| I I I I I I ■ I I I I I I I I I I I |_ 



1 2 

(H - K)vgg^ 



Fig. 5. — J — H vs. H — K color-color diagram (on the conventional, Vega-normalized system) 
comparing HDFN-JDl to a variety of cool or reddened stars. The NICMOS photometry has been 
converted to standard infrared magnitudes using synthetic color corrections, and is plotted as a 2a 
lower limit in J — H. The dotted lines mark the colors of a flat fiy spectrum (i.e., zero AB colors). 
The stehar data are taken from Frogel, Mould & Blanco 1990 and Costa & Frogel 1996 (LMC AGB 
stars), Whitelock et al. 1994, 1995 (IRAS Miras and carbon stars), Oestreicher, Schmidt-Kaler, 
& Wargau 1997 (supergiants), Frogel et al. 1990 (M giants), Leggett 1992 and Leggett et al. 1998 
(M dwarfs), Kirkpatrick et al. 1999 (L dwarfs), and Burgasser et al. 1999, Leggett et al. 1999, and 
Strauss et al. 1999 (methane dwarfs). IRC+10216 (Becklin et al. 1969) is also shown, along with 
LHS 3250 (Harris et al. 1999) which is believed to be a very cool DC white dwarf. The solid track 
shows the colors of a black body with labeled temperatures (in K). 
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Fig. 6. — Spectral energy distribution (SED) of HDFN-JDl. Open squares mark measurements 
with S/N < 2. Error bars (±lo") are superimposed, except when S/N < 1, where the downward 
error bar is shown as an arrow. Triangles mark 2a magnitude limits for C/300 through Jno- Three 
SEDs of other, real galaxies, artificially redshifted and normalized to the Kg flux of HDFN-JDl, are 
superimposed for comparison: the dusty galaxy HRIO shifted to 2; = 2.3, the HDF giant elliptical 
galaxy 4-752.1 shifted to 2: = 3.4, and the HDF Lyman break galaxy 2-901 shifted to z = 12.5. 
For the latter, the flux point shortward of Lyman a has been reduced to account for the assumed 
opacity increase of the Lyman a forest at z = 12.5, and flux shortward of the Lyman limit has been 
set to zero. For HDF 4-752.1, the shortest wavelength measurement (from the HDF C/300 image) 
probably overestimates the real UV flux from the gE galaxy because of a red leak in the WFPC2 
F300W filter; this is indicated by a downward-pointing arrow. All three SEDs are qualitatively 
similar to that of HDFN-JDl. 
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Fig. 7. — Infrared colors (on the AB system) versus redshift for galaxies with spectroscopic redshifts 
in the HDF-N. Point sizes scale with the -ffieo magnitudes. The solid horizontal lines mark the 
nominal colors of HDFN-JDl, and the dotted lines mark the la uncertainties for Hiqq — K and 
lower color limit for Jhq — Hiqq; the dashed line for Jno — -ffieo marks the 2a lower color limit. 
The curves show the expected colors of various model galaxy types. The reddest of these is a 
Bruzual & Chariot (1996) model for an elliptical galaxy formed at z = 15 in a 10^ year burst and 
evolving passively in a cosmology where Qm = 0.3, Qa = 0.7, and Hq = 70 km s~^ Mpc~^. The 
other curves are unevolving models matching the colors of local spirals and unreddened starburst 
galaxies. The evolving elliptical model matches the -ffieo ~ ^ color of HDFN-JDl at 3 ^ 2; < 6, 
but never quite reaches the 2a limit on Jno — Hiqq. The HI opacity of the intergalactic medium has 
been included by extrapolating the models of Madau (1995) to z = 15. At these large redshifts, the 
IGM transmission is effectively a step function at Lyman a. Redshifts were compiled from Cohen 
et al. (1996), Steidel et al. (1996b), Lowenthal et al. (1997), Dickinson (1998), Hogg et al. (1998), 
Weymann et al. (1998) and Spinrad et al. (1998), and Barger et al. (1999). 



